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From the observations of the anisotropies of the Cosmic Microwave Background 
(CMB) radiation, the WMAP satellite has provided a determination of the bary- 
onic density of the Universe, Q(,/i 2 , with an unprecedented precision. This imposes 
a careful reanalysis of the standard Big— Bang Nucleosynthesis (SBBN) calcula- 
tions. We have updated our previous calculations using thermonuclear reaction 
rates provided by a new analysis of experimental nuclear data constrained by R- 
matrix theory. Combining these BBN results with the Q.^h 2 value from WMAP, 
we deduce the light element ( 4 He, D, s He and 7 Li) primordial abundances and 
compare them with spectroscopic observations. There is a very good agreement 
with deuterium observed in cosmological clouds, which strengthens the confidence 
on the estimated baryonic density of the Universe. However, there is an important 
discrepancy between the deduced 7 Li abundance and the one observed in halo stars 
of our Galaxy, supposed, until now, to represent the primordial abundance of this 
isotope. The origin of this discrepancy, observational, nuclear or more fundamental 
remains to be clarified. The possible role of the up to now neglected 7 Be(d,p)2a 
and 7 Be(d,o) 5 Li reactions is considered. 



1. Introduction 

Big-Bang nucleosynthesis used to be the only method to determine the 
baryonic content of the Universe. However, recently other methods have 
emerged. In particular the analysis of the anisotropies of the cosmic mi- 
crowave background radiation has provided fi^/i 2 values with ever increasing 
precision. (As usual, Sib is the ratio of the baryonic density over the critical 
density and h the Hubble constant in units of 100 km-s _1 -Mpc _1 .) The 
baryonic density provided by WMAP 1 , fl b h 2 = 0.0224±0.0009, has indeed 
dramatically increased the precision on this crucial cosmological parameter 
with respect to earlier experiments: BOOMERANG, CBI, DASI, MAX- 
IMA, VSA and ARCHEOPS. It is thus important to improve the precision 
on SBBN calculations. Within the standard model of BBN, the only re- 
maining free parameter is the baryon over photon ratio 77 directly related to 
Q,bh 2 [£!b/i 2 =3.6519x 10 7 rj\. Accordingly, the main source of uncertainties 
comes from the nuclear reaction rates. In this paper we use the results of 
a new analysis 2,3 of nuclear data providing improved reaction rates which 
reduces those uncertainties. 



2. Nuclear reaction rates 

In a previous paper 4 we already used a Monte-Carlo technique, to calculate 
the uncertainties on the light element yields ( 4 He, D, 3 He and 7 Li) related 
to nuclear reactions. The results were compared to observations that are 
thought to be representative of the corresponding primordial abundances. 
We used reaction rates from the NACRE compilation of charged particles 
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reaction rates 5 completed by other sources 6 ' 7,8 as NACRE did not include 
all of the 12 important reactions of SBBN. One of the main innovative fea- 
tures of NACRE with respect to former compilations 9 is that uncertainties 
are analyzed in detail and realistic lower and upper bounds for the rates 
are provided. However, since it is a general compilation for multiple appli- 
cations, coping with a broad range of nuclear configurations, these bounds 
had not always been evaluated through a rigorous statistical methodology. 
Hence, we assumed a simple uniform distribution between these bounds for 
the Monte-Carlo calculations. Other works 10,11 have given better defined 
statistical limits for the reaction rates of interest for SBBN. In these works, 
the astrophysical S'-factors (see definition in Rcf. 4 ) were either fitted with 
spline functions 10 or with NACRE S-factor fits and data but using a dif- 
ferent normalization 11 . In this work, we use a new compilation 2 specifically 
dedicated to SBBN reaction rates using for the first time in this context 
nuclear theory to constrain the ^-factor energy dependences and provide 
statistical limits. The goal of the i?-matrix method 12 is to parametrize some 
experimentally known quantities, such as cross sections or phase shifts, with 
a small number of parameters, which are then used to interpolate the cross 
section within astrophysical energies. The i?-matrix theory has been used 
for many decades in the nuclear physics community (see e.g. Ref. 13 > 14 for 
a recent application to a nuclear astrophysics problem) but this is the first 
time that it is applied to SBBN reactions. This method can be used for 
both resonant and non-resonant contributions to the cross section. (See 
Rcf. 2 and references therein for details of the method.) The i?-matrix 
framework assumes that the space is divided into two regions: the internal 
region (with radius a), where nuclear forces are important, and the ex- 
ternal region, where the interaction between the nuclei is governed by the 
Coulomb force only. The physics of the internal region is parameterized by 
a number N of poles, which are characterized by energy E\ and reduced 
width j\. Improvements of current work on Big Bang nucleosynthesis es- 
sentially concerns a more precise evaluation of uncertainties on the reaction 
rates. Here, we address this problem by using standard statistical methods 
15 . This represents a significant improvement with respect to NACRE 5 , 
where uncertainties are evaluated with a simple prescription. The i?-matrix 
approach depends on a number of parameters, some of them being fitted, 
whereas others are constrained by well determined data, such as energies 
or widths of resonances. As usual, the adopted parameter set is obtained 
from the minimal x 2 value. The uncertainties on the parameters are evalu- 
ated as explained in Ref. 15 . The range of acceptable pi values is such that 
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X 2 (Pi) < X 2 (pT m ) + ^x 2 , where p™ m is the optimal parameter set. In this 
equation, A% 2 is obtained from P{y/2, A% 2 /2) = 1— p, where v is the num- 
ber of free parameters pi, P{a, x) is the Incomplete Gamma function, and p 
is the confidence limit (p = 0.683 for the la confidence level) 15 . This range 
is scanned for all parameters, and the limits on the cross sections are then 
estimated at each energy. As it is well known, several reactions involved in 
nuclear astrophysics present different data sets which are not compatible 
with each other. An example is the 3 He(a, -f) 7 Be reaction where data with 
different normalizations are available. In such special procedure 

is used 2 . 

This new compilation 2 provides \-a statistical limits for each of the 
10 rates: 2 H(p, 7 ) 3 He, 2 H(d,n) 3 He, 2 H(d,p) 3 H, 3 H(d,n) 4 He 3 H(a, 7 ) 7 Li, 
3 Hc(n,p) 3 H, 3 Hc(d,p) 4 He, 3 He(a, 7 ) 7 Be, 7 Li(p,a) 4 He and 7 Bc(n,p) 7 Li. The 
two remaining reactions of importance, n^->p and 1 H(n,7) 2 H come from the- 
ory and are unchanged with respect to our previous work 4 . 

3. SBBN calculations 

We performed Monte-Carlo calculations using Gaussian distributions with 
parameters provided by the new compilation and calculated the i He 1 D, 
3 He and 7 Li yield range as a function of r?, fully consistent with our pre- 
vious analysis 4 . The differences with Ref. 11 on the 7 Li yield is probably 
due to their different normalization procedure of the NACRE 5-factors. 
Figure 1 displays the resulting abundance limits (1-ct) [it was 2-a in Fig. 4 
of Ref. 4 ] from SBBN calculations compared to primordial ones inferred 
from observations. Using these results and the WMAP f^/i 2 range (quoted 
WMAP+SBBN in the following), it is now possible to infer the primordial 
4 He, D, 3 He and 7 Li abundances. 

We obtain (WMAP+SBBN) a deuterium primordial abundance of D/H 
= (2.601q'}7) x 10~ 5 [ratio of D and H abundances by number of atoms] 
which is in perfect agreement with the average value {2.78t°oii) x 10 ~ 5 01 
D/H observations in cosmological clouds 16 . These clouds at high rcdshift 
on the line of sight of distant quasars are expected to be representative 
of primordial D abundances. The exact convergence between these two 
independent methods is claimed to reinforce the confidence in the deduced 
flbh 2 value. 
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Figure 1. Abundances of 4 He (mass fraction), D, 3 He and 7 Li (by number relative to 
H) as a function of the baryon over photon ratio r\ or fl^h 2 . Limits (1-cr) are obtained 
from Monte Carlo calculations. Horizontal lines represent primordial 4 He, D and 7 Li 
abundances deduced from observational data (see text). The vertical stripes represent the 
(68% c.l.) Qf,h 2 limits provided by WMAP 1 or deduced from 7 Li and 4 He observations 
and SBBN calculations. For the dash-dotted lines in the bottom panel: see text. 
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The other WMAP+SBBN deduced primordial abundances are Yp = 
0.2457±0.0004 for the 4 He mass fraction, 3 He/H = (1.04 ± 0.04) x 10~ 5 
and 7 Li/H = (4.15^;4g) x 10~ 10 . We leave aside 3 He whose primordial 
abundance cannot be reliably determined because of its uncertain rate of 
stellar production and destruction 17 . 

The 4 He primordial abundance, Y p (mass fraction), is derived from 
observations of metal-poor extragalactic, ionized hydrogen (H II) regions. 
Recent evaluations gave a relatively narrow ranges of abundances: Y p = 
0.2452±0.0015 (Izotov et al. 18 ), 0.2391±0.0020 (Luridiana et al. 19 ). How- 
ever, recent observations by Izotov and Thuan 20 on a large sample of 82 
H II regions in 76 blue compact galaxies have lead to the value of Y p = 
0.2421±0.0021 that we adopt here. With this range, WMAP and SBBN 
results are hardly compatible. Nevertheless, as systematic uncertainties 
may prevail due to observational difficulties and complex physics 21 4 He 
alone is unsufficient to draw a conclusion. 

The 7 Li abundance measured in halo stars of the Galaxy is considered 
up to now as representative of the primordial abundance as it display a 
plateau 32 as a function of metallicity (sec definition in Ref. 4 ). Recent 
observations 31 have lead to (95% c.l.) Li/H = (1.23±g;|l) x 10~ 10 . These 
authors have extensively studied and quantified the various sources of un- 
certainty : extrapolation, stellar depletion and stellar atmosphere parame- 
ters. This Li/H value, based on a much larger number of observations than 
the D/H one was considered 4 as the most reliable constraint on SBBN and 
hence on Qbh 2 . However, it is a factor of 3.4 lower than the WMAP+SBBN 
value. Even when considering the corresponding uncertainties, the two 
Li/H values differ drastistically. This confirms our 4 and other 11,22 previ- 
ous conclusions that the Sl^h 2 range deduced from SBBN of 7 Li are only 
marginally compatible with those from the CMB observations available by 
this time (BOOMERANG, CBI, DASI and MAXIMA experiments). It is 
strange that the major discrepancy affects 7 Li since it could a priori lead 
to a more reliable primordial value than deuterium, because of much higher 
observational statistics and an easier extrapolation to primordial values. 

Fig. 2 shows a comparison between flth 2 ranges deduced either from 
SBBN or WMAP. The curves represent likelihood functions obtained from 
our SBBN calculations and observed deuterium 16 , helium 20 and lithium 31 
primordial abundances. These were obtained as in our previous analysis 4 
except for the new reaction rates and new D and 4 He primordial abun- 
dances. The incompatibility between the D and 7 Li likelihood curves is 
more obvious than before due to the lower D/H adopted value (Kirkman 
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Figure 2. Likelihood functions for D, 4 He and 7 Li (solid lines) obtained from our SBBN 
calculations and Kirkman et al. 16 , Izotov and Thuan 20 , and Ryan et al. 31 data for D, 
4 He and 7 Li respectively. The dashed curve represent the likelihood function for 4 He 
and 7 Li while the verical stripe shows the WMAP £2j,/i 2 range 1 . 

et al., averaged value). On the contrary, the new 4 He adopted value 20 is 
perfectly compatible with the 7 Li one as shown on Fig. 2 (likelihood curves) 
and Fig. 1 (abundances). Putting aside, for a moment, the CMB results 
on the baryonic density, we would deduce the following 68% c.l. inter- 
vals: 1.85< 7710 <3.90 [0.007<£Vi 2 <0.014] from 7 Li only or 5.4< 7710 <6.6 
[0.020<r2&/7 2 <0.024] from D only. If we now consider 4 He together with 
7 Li, we obtain 2.3< 7710 <3.5 [0.009<f2&/7 2 <0.013]. Hence, including these 
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new 4 He observations favors a low fl^h 2 interval as proposed in our previous 
work . The WMAP result on the contrary definitively favors the upper (D) 
one. If we now assume that the A He constrain is not so tight, because e.g. 
of systematic errors on this isotope whose weak sensitivity to Vlbh 2 requires 
high precision abundance determinations, the origin of the discrepancy on 
7 Li remains a challenging issue very well worth further investigations. 

4. Possible origins of 7 Li discrepancy between SBBN and 
CMB 

4.1. Stellar 

Both observers and experts in stellar atmospheres agree to consider that 
the abundance determination in halo stars, and more particularly that of 
lithium requires a sophisticated analysis. The derivation of the lithium 
abundance in halo stars with the high precision needed requires a fine 
knowledge of the physics of stellar atmosphere (effective temperature scale, 
population of different ionization states, non LTE (Local Thermodynamic 
Equilibrium) effects and 1D/3D model atmospheres 35 . However, the 3D, 
NLTE abundances are very similar to the ID, LTE results, but, neverthe- 
less, 3D models are now compulsory to extract lithium abundance from 
poor metal halo stars 36 . 

Modification of the surface abundance of Li by nuclear burning all along 
the stellar evolution has been discussed for a long time in the literature. 
There is no lack of phenomena to disturb the Li abundance: rotational 
induced mixing, mass loss,... 37 ' 38 . However, the flatness of the plateau over 
three decades in metallicity and the relatively small dispersion of data rep- 
resent a real challenge to stellar modeling. In addition, recent observations 
of 6 Li in halo stars (an even more fragile isotope than 7 Li) constrain more 
severely the potential destruction of lithium 39 . 

4.2. Nuclear 

Large systematic errors on the 12 main nuclear cross sections are 
excluded 2,3 . However, besides the 12 reactions classically considered in 
SBBN, first of all the influence of all nuclear reactions needs to be 
evaluated 3 . It is well known that the valley shaped curve representing 
Li/H as a function of rj is due to two modes of 7 Li production. One, at 
low r\ produces 7 Li directly via 3 H(a, 7) 7 Li while 7 Li destruction comes 
from 7 Li(p,a) 4 He. The other one, at high 77, leads to the formation of 7 Be 
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through 3 He(a,7) 7 Be while 7 Be destruction by 7 Be(n,p) 7 Li is inefficient 
because of the lower neutron abundance at high density ( 7 Be later decays 
to 7 Li). Since the WMAP results point toward the high rj region, a pecu- 
liar attention should be paid to 7 Be synthesis. In particular, the 7 Bc+d 
reactions could be an alternative to 7 Be(n,p) 7 Li for the destruction of 7 Be, 
by compensating the scarcity of neutrons at high r\. Fig. 1 shows (dash- 
dotted lines) that an increase of the 7 Be(d,p)2 4 He reaction rate by factors 
of 100 to 300 would remove the discrepancy. The rate for this reaction 9 can 
be traced to an estimate by Parker 40 who assumed for the astrophysical 
S-factor a constant value of 10 5 kev.barn. based on the single experimen- 
tal data available 41 . To derive this S'-factor, Parker used this measured 
differential cross section at 90° and assumed isotropy of the cross section. 
Since Kavanagh measured only the po and pi protons (i.e. feeding the 8 Bc 
ground and first excited levels), Parker introduced an additional but ar- 
bitrary factor of 3 to take into account the possible population of higher 
lying levels. Indeed, a level at 11.35 MeV is also reported 42 . This factor 
should also include the contribution of another open channel in 7 Bc+d: 
7 Be(d,a) 5 Li for which no data exist. In addition, one should note that no 
experimental data for this reaction is available at energies relevant to Be 
Big Bang nucleosynthesis (Fig. 3), taking place when the temperature has 
dropped below 10 9 K. A seducing possibility 3 to reconciliate, SBBN, 7 Li 
and CMB observations would then be that new experimental data below 
E d = 700 keV (E cm w0.5 MeV) for 7 Bc(d,p)2 4 He [and 7 Be(d,a) 5 Li] would 
lead to a sudden increase in the S'-factor as in 10 B(p,a) 7 Be 43 ' 5 . This is 
not supported by known data, but considering the cosmological or astro- 
physical consequences, this is definitely an issue to be investigated and an 
experiment is planned in 2004 at the Cyclotron Research Centre in Louvain- 
la-Neuve. 

4.3. Cosmology 

Recent theories that could affect BBN include the variation of the fine struc- 
ture constant 44 , the modification of the expansion rate during BBN induced 
by quintessence 45 , modified gravity 46 , or leptons asymmetry 47 . However, 
their effect is in general more significant on A He than on 7 Li. 

It may not be excluded that some bias exists in the analysis of CMB 
anisotropies. For instance, it has been argued 48 that a contamination of 
CMB map by blazars could affect the second peak of the power spectrum 
on which the CMB flbh 2 values are based. 
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Figure 3. The only experimental data available for the 7 Be(d,p)2 4 H reaction from Ka- 
vanagh (1960). The displayed S-factor is calculated as in Parker (1972) from the differ- 
ential cross section at 90° (x47r) leading to the ground and first 8 Be excited states. Note 
that no data is available at SBBN energies as shown by the Gamow peak for a typical 
temperature of Tg = 0.8 



4.4. Pregalactic evolution 

We note that between the BBN epoch and the birth of the now observed 
halo stars, sal Gyr have passed. Primordial abundances could have been 
altered during this period. For instance, cosmological cosmic rays, assumed 
to have been born in a burst at some high redshift, could have modified 
these primordial abundances in the intergalactic medium 49 . This would 
increase the primordial 7 Li and D abundances trough spallative reactions, 
increasing in the same time the discrepancy between SBBN calculations 
and observations instead to reconciliatc them. 

Another source of alteration of the primordial abundances could be the 
contribution of the first generation stars (Population III). However, it seems 
difficult that they could reduce the 7 Li abundance without affecting the D 
one, consistent with CMB fl^h 2 . 
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5. Conclusions 

The baryonic density of the Universe as determined by the analysis of the 
Cosmic Microwave Background anisotropies is in very good agreement with 
Standard Big-Bang Nucleosynthesis compared to D primordial abundance 
deduced from cosmological cloud observations. However, it strongly dis- 
agrees with lithium observations in halo stars (Spite plateau) and possibly 
4 He new observations. The origin of this discrepancy if not nuclear, is a 
challenging issue. 
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